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Sound generation and -interaction is highly complex, nonlinear and self-organized. Already 150 
years ago Lord Rayleigh raised the following problem: Two nearby organ pipes of different fundamen- 
tal frequencies sound together almost inaudibly with identical pitch. This effect is now understood 
qualitatively by modern synchronization theory (M. Abel et al., J. Acoust. Soc. Am., 119(4), 2006). 
For a detailed, quantitative investigation, we substituted one pipe by an electric speaker. We ob- 
serve that even minute driving signals force the pipe to synchronization, thus yielding three decades 
of synchronization - the largest range ever measured to our knowledge. Furthermore, a mutual 
silencing of the pipe is found, which can be explained by self-organized oscillations, of use for novel 
methods of noise abatement. Finally, we develop a specific nonlinear reconstruction method which 
yields a perfect quantitative match of experiment and theory. 
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Introduction. In a seminal publication Lord Rayleigh 
reported on experiments involving two organ pipes of 
close pitch, positioned next to each other. He observed 
the following peculiar behavior: alone, each pipe sounded 
with its own natural frequency. Together they sounded 
in perfect unison and almost reduced one another to si- 
lence [T]. This phenomenon can be described by synchro- 
nization theory, a general nonlinear principle with strik- 
ing applications in the natural sciences including physics, 
chemistry, neurology or biology |2]. 

Here, we focus on an organ pipe as i) prototypical 
for aero-acoustical sound generation, ii) paradigmatic for 
synchronization of coupled oscillators and iii) a beautiful 
musical instrument for which we present a mathemati- 
cal model. By driving one organ pipe sinusoidally by an 
electric speaker we can demonstrate in great detail the 
importance of nonlinear effects in sound generation and 
-interaction [3 in contrast to linear response theory [4^ 
conventionally applied for musical instruments. 

The general description of the dynamics of an organ 
pipe is given by the compressible Navier-Stokes Equa- 
tions with suitable boundary conditions. One can solve 
the equations numerically , or investigate them analyt- 
ically [6^ . Both ways reproduce different aspects of sound 
production. Here, we are interested in the interaction of 
an organ pipe with a sound source. Then, it is of advan- 
tage to model only the relevant characteristics in terms of 
reduced models [7J. Such an elementary model allowing 
for complex dynamics is given by an autonomous oscil- 
lator [2i i8j, which includes an oscillatory unit, energy 
supply, and energy loss by radiation and damping. 

Let us identify these units in the organ pipe. Energy 
is supplied steadily by the wind system through the pipe 
foot and establishes a turbulent vortex street. Each time 
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a vortex detaches, a pressure fluctuation enters the res- 
onator, inside which characteristic waves are selected, 
and radiated at the pipe mouth by an oscillating air- 
sheet [3 [3] . In our model, this air-sheet constitutes the 
basic oscillating unit. Inside the resonator Sound Pres- 
sure Levels (SPLs) up to 160 dB can occur, such that 
viscous damping contributes in energy dissipation. Ex- 
ternal acoustical fields couple to the system through the 
air-sheet, possibly described by a (nonlinear) acoustical 
admittance [lOj. By Lighthills analogy [6 , a coupling 
by the turbulent vortex street is expected to be of lower 
order. 

The above scenario can be described by a reduced, two- 
dimensional model for the oscillatory unit ^: 

l-5(e,o-o, (1) 

where the function §(^,0 contains the above mentioned 
ingredients and with the condition that a limit-cycle solu- 
tion of frequency exists. Then, with ^ ^(i)-e"^*^*^ the 
phase (j) and the amplitude A are well defined. An exter- 
nal driving enters on the right hand side of Eq. ([T|) . Cor- 
responding to the experiment, sinusoidal driving is used: 
e sin(27rz/t -I- 0o) with e the coupling strength. Here, one 
assumes that an oscillator represents the basic physics of 
the pipe with regard to sound generation and synchro- 
nization - typically, the oscillating air sheet [7]. Because 
the air sheet is the source of sound radiation, the mea- 
surement at the microphone can be taken as the state of 
the oscillator (with a phase shift accounting for distance). 

Close to the limit cycle the amplitude is slaved by the 
phase, allowing the description in terms of the phase dif- 
ference, W, between driving and oscillator 

<F = ~2TT{i^~i'o) + eq{<F). (2) 

The parameters are driving frequency v, and coupling e. 
The study of the parameter plane (v, e) yields triangular- 
shaped synchronization regions, the well-known Arnold- 
Tongues im. 
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FIG. 1: Sketch of the experiment. Pipe and loudspeaker stand 
side by side, the signal is measured by a microphone with 
equal distance to both sound sources. 



In order to determine Eq. ([T]) directly from data, we 
have elaborated a numerical method based on embed- 
ding theory \12\ which allows for a reconstruction and 
comparison of the characteristics of model and data. We 
recover the power spectra and synchronization properties 
of the organ as an acoustical systems. This opens the way 
for a detailed theoretical investigation of the system and 
hints to how a model can be derived from first principles. 

The experiment. The setup of the experiment is 
sketched in Fig. 1. The pipe |T3] was wooden and closed 
at the upper end, tuned at = 686Hz. It was driven 
by an especially fabricated miniature organ 13J with a 
blower connected to the wind-belt and further by flex- 
ible tubes to the wind-chest. Measurements took place 
inside a suitable anechoic box. The loudspeaker for the 
external, driving sound signal, was positioned side-by- 
side to the organ-pipe, cf. Fig. 1. The emitted signals 
were registered at a distance of 16 mm to either pipe and 
speaker, allowing for constructive or destructive interfer- 
ence of the superposed signals. To ensure that the phase 
of the pipe is correctly detected, we carried out simulta- 
neous control measurements inside and outside the pipe 
and at the microphone - all results were consistent. 

To explore the coupling-detuning plane {v, e), the loud- 
speaker SPL and frequency were varied separately; the 
first between -l-lOdB and -50 dB relative to the reference 
signal of the organ pipe in steps of 2 dB, the latter accord- 
ing to the size of the synchronization range. To determine 
the minimal achievable resolution in z/, we consider the 
sources of variations of the pipe's frequency. The wind 
pressure was 700±9 Pa, giving a frequency variation of 
±0.1 Hz at the scale of seconds; the temperature varied 
with the circadian rhythm at 292±1 K with a resulting 
variance of ~2 Hz [14] . This slow change did not affect an 
individual run with fixed SPL., however, runs with dif- 
ferent SPL were adjusted to the circadian variation. We 
measured the synchronization range down to the maxi- 
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FIG. 2: Synchronization plot: if the detuning of pipe and 
loudspeaker is small enough, frequency synchronization is 
found as a plateau. Within the synchronization region, the 
phase difference iP" varies over an interval of size tt, where 
the measured SPL shows a sharp minimum for 'I' — n. Top 
[SPLspeaker ~ SPLpipe): the frequency shows a very nice 
plateau with a saddle-node bifurcation; for 'I' — n nega- 
tive interference is observed by a decrease of 6 dB. Bottom 
{SPLspeaker — SPLpipe -\- 10 dB): the transition from syn- 
chronization is quite sharp, indicating strong nonlinearities. 
The amplitude decrease can be explained by synchronization 
with an additional Helmholtz-resonator-like behavior of the 
pipe, deepening the SPL gap to 15 dB. The maximum is ob- 
tained by addition of the amplitudes to 12 dB. 



mal resolution of 0.1 Hz as set by the wind supply; below 
noise destroys synchronization. With respect to noise, we 
enhanced the signal-to-noise ratio was enhanced by long- 
time averaging, such that irregular phase slips are leveled 
out and a synchronization region can be obtained even 
for very small driving. This was achievable by heavily 
automatized measurements of a total duration of about 
3 weeks. 

We investigated two acoustically relevant character- 
istics: the frequency difference, Ajz/i/q = {vspeaker — 
vpipe) /t'oi and the spectrum of the measured loudspeaker 
signal. Their dependence on the detuning is shown in 
Fig. 2 for two exemplary relative driving strengths, dB 
and 10 dB. Note, that vpipe is the frequency of the driven 
pipe which might be different from its natural frequency. 

The transition to synchronization can be seen from the 
graph for the frequency difference. For equal SPL and 
below (upper panel in Fig. 2), a saddle node bifurcation 
is found, as predicted theoretically [SJ [TS] . For couplings 
stronger than OdB, the bifurcation tends to become very 
sharp, indicating that the weakly nonlinear approxima- 
tion [2j breaks down (lower panel Fig. 2) . 

In the graphs for the amplitude, one recognizes two 
effects: synchronization and resonance. Due to synchro- 
nization the phase shift between the two emitted signals 
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FIG. 3; Arnold tongue measured over approximately 3 
decades. The left plot shows the Tongue obtained from the ex- 
perimental data, while the right one shows the reconstructed 
tongue. The color coding corresponds to the amplitude mea- 
sured at the microphone position at the loudspeaker frequency 
Sifspeaker)- To guide the eye, the synchronization region is 
marked by a straight line. The synchronization edge grows 
linearly with the coupling, seen as a logarithmic bend in the 
used semi-log plot. The lowest possible frequency resolution, 
0.1 Hz, corresponds to the variation of the frequency produced 
by the wind supply, shown by the inset. The reconstructed 
Arnold tongue is much more symmetric. This indicates again 
strong nonlinearities. The coincidence between experiment 
and model is almost perfect in the low-coupling region. 



varies within an interval of size tt; as for Hclmholtz- 
resonance, the loudspeaker signal is phase-shifted by ~ tt 
and re-emitted by the pipe. From synchronization alone, 
the superposition of the signals of speaker and pipe can 
only vary within SPLspeaker i SPLpip,, ; with resonance 
much weaker amplitudes result, as seen in Fig. 2. This re- 
sult implies a novel way of sound reduction, where the re- 
ductor acts as an active element, adjusting the frequency 
without any external control. 

The full variation in parameter space is shown in the 
left panel of Fig. 3 by an Arnold Tongue, here plotted log- 
arithmically due to the enormous coupling range investi- 
gated. The color coding corresponds to the the spectral 
power of the loudspeaker frequency S^vspeaker), mea- 
sured at the microphone. We varied the amplitude of the 
loudspeaker from -1-10 dB to -50 dB, relative to the SPL 
of the organ pipe, in steps of 2dB. The range covers 3 
decades - the widest range ever measured in synchroniza- 
tion experiments. The synchronization edges are marked 
in Fig. 3 by white lines. Clearly, the linear shape is bent 
due to the semi- log plotting. 

Reconstruction of the dynamical system. We do not 
have access to the "state" of the oscillating air-sheet, 
but we know the recorded SPL, x{t), of the organ pipe. 
By embedding theory we can infer a differential embed- 
ding {x,x,x, . . .) with a maximal embedding dimension 
of three [12] |T5] . In this space there exists an equivalent 
to Eq. ([T|: X — f{x,x) = 0. We reconstruct this equa- 
tion step by step to further investigate its dynamical and 



predictive properties. 

The data series consists of 110250 data points, with a 
sampling rate of 11025 Hz. Normalization in space and 
time yields variance 0.5 and frequency 1. The crucial 
computation of derivatives was accomplished by spec- 
tral smoothing |17j : i) Fourier transformation, ii) 8th 
order Butterworth-filtering with cutoff at 4.5 i^o to sup- 
press noise amplification , and Hi) back-transformation. 
If the cutoff is increased more harmonics enter the filtered 
time series, if the cutoff is too low too much information 
is filtered out, i.e. the necessary nonlinearities are sup- 
pressed. Thorough testing yielded that in our case a three 
dimensional embedding do not improve the results. 

The unknown function / is estimated by nonparamet- 
ric regression, formulated as a minimization problem: 

\\x — f{x, i)\\2 — min, with || • ||2 the /^-norm of the data 
vector. The unknown function / is found by variation 
in function space, where we used, for the sake of com- 
putational simplicity, polynomials of order three [TH] HH] ; 
higher orders do not improve the model. Specifically, 
/(x, x) = ao + aix + 02^^ -t- a^x^ + a^x + a^x'^ + aex^ + 
a^xx + a^x'^x + agxx^ with aq — 0.25, ai — —0.92, 02 ~ 
-0.18,03 = -0.12, a4 = 0.20,05 = -0.33,06 = 
0.056,07 = -0.015,08 = -0.923,09 = -0.072 . Note 
the striking similarity to the van der Pol oscillator with 
fvdP — —X + x{\ — a;^), reflected in the dominant coeffi- 
cients oi and og. The latter is responsible for nonlinear 
damping, whereas energy is supplied by the constant and 
the i-term. Other terms assist nonlinear damping and 
are indispensable to find the correct frequencies in the 
reconstructed system. Since the observables are not di- 
rectly related to the physical driving (wind) and damping 
mechanisms we hesitate giving a complete physical inter- 
pretation. 

Numerical stability analysis yields a repelling fixed 
point at (0.254,0) and an attracting limit cycle, plot- 
ted in the inset of Fig. 4, which shows convincing coinci- 
dence with the filtered experimental data. For acoustical 
comparison, we compare the power spectra of pipe and 
reconstruction in Fig. 4. The positions of the harmonics 
are in perfect agreement, and their ratio coincides well. 

Finally, the model is synchronized with an external, 
sinusoidal driving. For the coupled equation, we solve 
X — f{x, x) + esin(cLii), with the driving frequency lu and 
the coupling parameter e = 0.025 • SPLspeaker- The 
latter relation describes the infiuence of the loudspeaker 
on the pipe with a factor 0.025, determined in order to 
obtain optimal coincidence of experimental and recon- 
structed Arnold Tongue, see Fig. 3. Recently, a way to 
extract the coupling from data has been proposed |20j : 
in principle it can be obtained by a detailed analysis of 
acoustics and fiuid dynamics at the pipe mouth. 

Conclusion. Since the time of Lord Rayleigh the non- 
linear interaction of acoustical sources is under discussion 
[Ij. We highlight the acoustical effects of phase synchro- 
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FIG. 4: Left: The power spectra of the measured (black) 
and the reconstructed signal (blue). The musical sound is 
reproduced quite nicely, as recognized by the coincidence of 
the spectra; for better visibility the maxima of the measured 
signal are shown by black dots. The inset shows the time 
series and its reconstruction in the embedding space, which 
coincide very well. Right: Contour plot of the function / 
found by nonparametric regression, the color coding is shown 
by the color bar on the right. 



nization, experimentally realized by an organ pipe driven 
externally by a loudspeaker. The application of synchro- 
nization theory suggests a novel type of sound control, 
where the passive element adjusts its frequency exactly to 
the source. Since the system involved can be abstracted, 
we conclude that such a control can be applied to a di- 
versity of situations, from musical instruments to noise 
reduction in vibrating systems, be it mechanical or hy- 
dromechanical. 

With respect to synchronization, we found the deepest 
Arnold tongue ever seen experimentally, suggesting wind 
instruments as paradigmatic for synchronization. The re- 
sults confirm the theory for small coupling; for large cou- 
pling we give experimental access to nonlinear correction 
terms to be analyzed further. To analyze the acoustical 
properties of the pipe we propose an autonomous oscil- 
lator, reconstructed from a novel type of data analysis. 
The agreement between model and experiment in terms 
of Arnold Tongues and power spectra is excellent - mu- 
sical and synchronization characteristics are well repro- 
duced. Conventional methods, as transfer functions, or 
admittance [2T] do not allow such a direct interpretation 

In this work acoustics is paired with nonlinear dynam- 
ics and data mining. Organ builders have developed com- 
plicated empirical rules to arrange organ pipes within 
a register to avoid synchronization effects. Our results 
do not only allow for an easy simulation of instruments, 
but as well development cycles and tuning of instruments 
could be enhanced. Noise reduction is possible by self- 
organization of two sound sources such that they interfere 
negatively - no external control is needed. The applica- 
tions of this principle might be interesting for a variety 
of technically important situations. 
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